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Abstract 

General Relativity simplifies dramatically in the limit that the number of spacetime di- 
mensions D is infinite: it reduces to a theory of non-interacting particles, of finite radius 
but vanishingly small cross sections, which do not emit nor absorb radiation of any finite 
frequency. In many respects, black holes and black branes behave like configurations of 
dust. The simplicity of this limit motivates the study of the theory in an expansion inl/D. 
Large D introduces a separation of scales that allows an effective theory of black hole dy- 
namics. Wc develop to leading order in 1/D this effective description for massless scalar 
fields and compute analytically the scalar absorption probability. We solve to next-to- 
next-to-leading order the black brane instability, with very accurate results that improve 
on previous approximations with other methods. These examples demonstrate that prob- 
lems that can be formulated in an arbitrary number of dimensions may be tractable in 
analytic form, and very efficiently so, in the large D expansion. 
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1 Introduction 



The fascination power of General Relativity stems largely from the wealth of physical 
phenomena that are encoded in equations as simple as 



Naturally, this conceptual simplicity unfolds its rich dynamics at the cost of technical 
complexity. It is very difficult to find closed exact solutions to these coupled, non linear, 
partial differential equations, for almost any phenomenon of interest unless a substantial 
degree of symmetry is present. 

Physical theories often contain parameters that can be varied in such a way that the 
theories remain well defined. A fruitful strategy is to focus on regions of the parameter 
space, usually close to its origin or boundaries, where the theories simplify. Einstein's 
theory in vacuum, (1.1), appears to have only one natural parameter: the number D 
of spacetime dimensions. We will argue that in the limit D — t- oo. General Relativity 
simplifies dramatically, its dynamics becoming trivial at all non-zero length scales. This 
is a strong motivation for the study of the theory in an expansion in 1/D. 

While it seems unlikely that our universe be infinite-dimensional, the study of General 
Relativity around this limit can be useful, both to gain a better understanding of the theory 
and as an approximation scheme for calculations in less unrealistic cases, say D = 4 oi 
D = 10, 11. We have learned in recent years that new features appear in the spectrum 
of black hole solutions as D grows beyond four [1]. Still, it may not be unreasonable to 
expect that some properties remain more robust as D is increased — to begin with, the 
theory (1.1) does not have black holes nor a dynamical graviton when D <3, and always 
has them for any D > 4. The question of to what extent an expansion in 1/D is a good 
qualitative guide to moderate D, and if so, how accurate it is, probably depends on the 
specific problem under consideration. In this article this concern will remain mostly in 
the background, and instead we will focus on understanding the main properties of the 
limit and on how to organize calculations in the 1 /D expansion. Nevertheless, one of our 
examples shows that these techniques can give very accurate results even at relatively low 
values of D. 

Early studies of gravity in the large D limit analyzed the quantum theory and the 
properties of its Feynman diagrams [2, 3]. The main motivation is the possible analogy with 
the large limit of SU{N) gauge theories — indeed, the local Lorentz group S0{D — 1, 1) 
is the gauge group of gravity. The large N limit of Yang-Mills theories is useful because, 
although the number of gluons grows unbounded, they arrange themselves into worldsheets 
of strings (propagating in more dimensions). One might hope that a similar miracle could 
occur also for large D gravity. However, D appears not only in the number of graviton 
polarizations but also, more troublingly, in phase space integrals. Actually, taking D 
to be large seems a bad idea for a quantum field theory, since the ultraviolet behavior 
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generically worsens. This can be alleviated by focusing on Kaluza-Klein truncations of 
the spectrum [4, 5], which retain the growing number of polarizations but make the short- 
distance behavior essentially four- dimensional, with the usual divergence problems in the 
gravitational sector. We shall not pursue any of these approaches. 

Instead, we study mostly the classical theory of eq. (1.1), which is well defined in any 
D. After all, many quantum properties of gravity are dominated by classical effects such 
as black hole formation. One speculation is that the illnesses of quantum gravity may 
be absent outside the horizons when D — )• oo and might remain under some control close 
to the horizon in the 1/D expansion. At any rate we will see that, even if its quantum 
version happens to be badly behaved, the large D limit of the classical theory is a useful 
one. 

The study of classical gravity in arbitrary D > 4 has gained momentum over the 
years, and it has often seemed natural to examine specific results in the limit D — )■ oo. 
However, very rarely has the large D expansion been pursued as a subject in its own 
right. A notable exception are refs. [6, 7], which are closest in spirit and techniques to 
our approach. Nevertheless, their context was restricted to a particular phenomenon (the 
Euclidean zero mode of the Schwarzschild solutions) and a bigger framework was not 
developed. Some general observations about black holes and black branes at large D are 
made in [8, 9] emphasizing slightly different features than here. Although there are many 
studies of gravitational phenomena in arbitrary D, those that make explicit reference to 
the large D limit are, as far as we are aware, relatively few, e.g., [10]-[21]. 

In summary: a systematic approach to the large D limit of classical General Relativ- 
ity has been lacking so far. We aim to provide some basic entries to the concepts and 
techniques of this subject. 

A main element of our approach is that, in contrast to the studies inspired by large N 
gauge theories, we will not focus on perturbations around a Minkowski background but 
rather on non-perturbative objects in the theory, namely its black holes. We argue that 
in the limit D — >• oo black holes are in many respects non-interacting particles. They do 
not attract each other and, even if their radius remains finite, their collision cross sections 
vanish. They reflect perfectly all radiation of any finite frequency. Some of these features 
suggest to think of black holes as 'dust particles'. Black branes do behave like if made of 
dust, with no tension to hold them together. 

The simplicity of this limit makes it a good starting point for a perturbative expansion 
in 1/D. We will see that in the limit D ^ oo the gravitational field vanishes outside the 
horizon at r = tq, while at large but finite D this field is strongly localized close to the 
horizon in the region 

r-ro<g. (1.2) 

Crucially, a new scale appears owing to the very steep gradient, ~ D/vq, of the gravi- 
tational potential near the horizon. This separation between scales r^/D ^ tq allows to 
develop an effective theory of black hole dynamics. Fields outside the black hole propagate 
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in an effectively flat spacetime, subject to certain boundary conditions very near the hori- 
zons, which replace the black holes. Technically, this is a problem of matched asymptotic 
expansions. Conceptually, it is an effective theory in which the degrees of freedom for 
the black hole are integrated out and replaced by boundary conditions on large-distance 
fields. What is peculiar to this effective theory is that the notion of short-distance degrees 
of freedom is the result of having large -D, instead of the more conventional idea of con- 
sidering wavelengths much larger than the horizon radius. This results in a much larger 
range of applicability of the effective theory. 

We begin in the next section motivating the 'small, non-interacting particle' picture of 
the limit — )• oo through an extensive study of known black hole solutions. In section 3 
we introduce the notion of the 'sphere of influence' of the black hole. In section 4 we 
discuss how and when classical gravitational radiation can be emitted through black hole 
interactions at D — t- oo. In section 5 we make some comments about quantum effects and 
Hawking radiation in black holes at large D. In section 6 we introduce the 1/D expansion 
in the study of the propagation of massless scalars in the black hole background. We 
solve the theory in the region near the horizon and find the effective boundary conditions 
for outside fields. We obtain a compact analytic expression for the scalar absorption 
probability that is valid over a very wide range of frequencies. In section 7 we solve the 
spectrum of unstable perturbations of black branes in closed analytic form to next-to- 
next-to-leading order at large D. The results are very accurate for all but the two lowest 
dimensions. We conclude in section 8. 

2 Large D limit of black holes 

We take the point of view that General Relativity in vacuum, the theory of (1.1), is 
essentially a theory of black holes that interact via the gravitational field between them, 
and which can emit and absorb gravitational waves through these interactions. We can 
also include as objects of study the singular plane wave solutions that appear in the infinite 
boost limit of black holes, but in general we exclude other nakedly singular solutions. 
We begin with the most basic solution, the Schwarzschild-Tangherlini spacetime 



with a horizon of radius tq [22]. When we take the limit D — t- cxd we have to decide how 
the magnitudes in the problem scale with D, and in particular which quantities are kept 
fixed, i.e., do not scale with D. An appropriate choice is one that results in a simple limit. 
We will keep the length ro fixed as — t- oo. 

Although we could set tq equal to one, the discussion is often clearer if we keep it 
explicit. 



fdt^ + f-^dr"^ + r'^dQ.D-2 





(2.2) 
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2.1 'Smallness' of the horizon and a hierarchy of scales 

Some of the properties of black holes at large D are not due to spacetime curvature but 
rather follow from elementary flat-space geometry (see appendix A). In particular, when 
D grows large, the area of the unit-radius sphere S^"^, 

27r(^-i)/2 

vanishes as 



= ^ ,D-i\ ' '^^•^) 



We may then say that these spheres becomes increasingly small at very large D, shrinking 
at a rate ~ D~^l'^ . In particular, if we shoot a projectile at this unit-radius sphere, 
the impact cross section vanishes when D ^ oo, even if the projectile will hit the target 
whenever the impact parameter is < 1.^ This geometric effect will be present in all 
calculations of total cross sections at D — )• oo. The strong semifactorial suppression 
~ D~^/'^ may sometimes hide effects that one is interested in, and in these cases it may 
be convenient to eliminate it by, e.g., considering ratios of appropriate magnitudes. 

There are other properties of the horizon, independent of this flat-space effect, which 
indicate that it is useful to think of the large D black hole as a very small object. The 
surface gravity on the horizon is 

L> - 3 , , 

K=^— , 2.5 

which diverges as D — >• oo. Thus a length-scale arises associated to /t, 

4 = (2-6) 

which is parametrically much smaller than the horizon radius tq. The intrinsic curvature 
gives essentially this same length. The Kretschmann scalar in the geometry (2.1) is 

^w-^ '°-'"° -^' (^)''°"' (2.7) 

so the characteristic curvature length of the horizon is 

{K{ro))-"' ^ § ~ 4 • (2-8) 

These short scales essentially reflect the large radial gradient of the gravitational potential 
near the horizon, 

drfl, ^ ^ • (2.9) 

Since the classical theory of (1.1) is scale- free, one usually regards ro as the only length 
scale present in the solution (2.1), and in this respect it would not make sense to say that 



0,0-3/ {D ^ 2), whose dominant large D behavior is 
inversely semifactorial like that oi Q,d~2- 
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the black hole is small (compared to what?). But when D is very large we have identified 
a much shorter scale characterizing important aspects of the geometry. 

Summarizing, there are two senses in which we can say that the black hole is effectively 
small. One is the vanishing area effect (2.4). The other is the appearance of the small 
length scale ro/D, which we will see controls much of the classical black hole physics. 
Along this article we will develop and refine these ideas in much more detail. 

2.2 Absence of interactions 

Outside the horizon the gravitational potential (ro/r)^~^ vanishes exponentially fast in D. 
When Z) — 7- cxD the lines of force are infinitely diluted in the infinite number of directions 
available, so there is no gravitational force outside the horizon. The horizon itself becomes 
a surface of infinite curvature. The spacetime of the black hole is then a flat geometry 
with a sphere cut off at the radius r = tq. The area of this sphere is vanishingly small. 

This implies that at — )• oo there is no force at all between two black holes. Multiple 
black hole solutions are obtained by simply cutting off spheres at different places in flat 
spacetime. Given this absence of interactions, and that these black holes present a van- 
ishing cross section for colliding with each other, we may regard this as a system of 'dust 
particles'. 

The absence of a gravitational force, or equivalently the flatness of the metric outside 
the horizon, is also reflected in the fact that, for the solution (2.1), we have that 

{D-2)nn-2 n-3 (2.IO) 

vanishes as D — )• 00 due to the factor r2z)_2. A caveat is now in order. The classical 
theory of (1.1) is in essence a purely geometric theory that only contains geometrical 
magnitudes, such as lengths and areas. Any other physical magnitudes, such as mass 
or angular momentum, which require the introduction of conversion constants such as 
G (absent from (1.1)), have only secondary meaning. Thus, G and M do not have any 
independent meaning in vacuum gravity. Nevertheless, with a view to coupling gravity to 
other matter systems, it may seem more appropriate to interpret (2.10) as saying that the 
gravitational coupling vanishes like 

G ^ D-^/^ ^0, (2.11) 

which suggests that the gravitational field becomes a free field at large D. 

All large D scales in vacuum gravity must have an interpretation of purely geometric 
origin. GM can be regarded as a measure of how the extrinsic curvature of surfaces of 
constant r differs from the extrinsic curvature of the same surfaces when embedded in 
Minkowski space. At large D, this difference vanishes and hence GM — )■ 0. We may 
consider introducing a gravitational mass-length scale 

iM ~ (GM)i/(^-3) ~ ^ , (2.12) 
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characterizing this effect. This length, however, can equally well be obtained from the 
horizon area, 

iA - Af""-^^ ~ ^ . (2.13) 



It derives from the geometric effect discussed in eq. (2.4). The conceptual prevalence of 
geometric magnitudes must be borne in mind whenever we choose, for physical illustration, 
to frame our discussion in terms of non- geometric, secondary quantities.^ 

The picture of black holes as non-interacting dust also agrees well with other properties. 
The Bekenstein-Hawking entropy of the Schwarzschild black hole behaves like 

S{M) ~ M . (2.14) 

In contrast to the situation at finite -D, where S oc with a > 1, the fact that S (x M 
means that there is no entropic gain in merging two black holes. Nor is there any entropic 
penalty in splitting a black hole in two (recall that the horizon becomes singular in the 
limit). This is a reflection of the absence of interactions noted above. 
Consider now a black p-brane, 

p 

ds"^ = -fpdt^ + + fp^dr^ + r^dnD^p^2 , (2.15) 

i=l 

This brane is characterized by an energy density e and a pressure P along its worldvolume 
such that 

P = . (2.17) 

When D — )• oo this pressure vanishes: the brane has a dust equation of state. 

The instability of the dust brane to fragmentation in this limit is easy to establish. At 
any finite D, neutral black branes suffer from a Gregory-Laflamme instability to growing 
inhomogeneities along the worldvolume, and are expected to eventually break up [25]. The 
threshold mode at the onset of this instability has been studied at large in [6, 7], with 
the result that perturbations with wavelength longer than 

Ae = ^ (1 + 0(1)-^)) (2.18) 



are unstable. Thus, when D ^ oo perturbations of arbitrarily short wavelength drive the 
break up of the brane. We will revisit this instability in more detail in section 7. 

For black p-branes we can regard p as a parameter that can scale with D in different 
manners. The case of p ~ 0{D^) has been discussed above. Another possibility is that 



■^Note that the quantum theory introduces the Planck length, which may be chosen to scale or not with 
D, see sec. 5. 

^The shrinking effect of angular spheres may be eliminated by considering, e.g., the ratio between initial 
and final entropies. This also eliminates the dependence on G and the Planck length. 
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p ~ D in such a way that n = D — p — 3 remains finite, i.e., we keep fixed the codimension 
of the brane instead of its worldvolume dimension. The gravitational field of these branes 
remains finite outside the horizon, and so does, too, their energy density and pressure. 
It seems appropriate to regard these black branes as belonging in a different sector of 
the theory than the infinitely localized ones we are considering so far, as they cannot 
be produced through processes involving a finite number of the latter. Their dynamics 
appears to be dominated by the degrees of freedom of a Kaluza-Klein reduction down to 
n + 3 dimensions. Except for a brief mention in the next subsection, we will not consider 
these branes in the remainder. 

2.3 Other black holes: rotation and other topologies 

When rotation is present in any D > 4 there are large classes of black hole solutions with 
many possible horizon topologies. Although the entire space of solutions has not been 
fully mapped, we can argue that the limit D — )• oo results in configurations with singular 
horizons and flat space outside them. Equilibrium horizons are equipotential surfaces, and 
at large D the gravitational potential falls off steeply in the transverse direction away from 
them. If we take the limit D — t- oo in such a way that the length scales that characterize 
the horizon remain flnite, then the potential will drop infinitely fast outside the horizon. 
A possible concern is that in D > 4 there are horizons that are much more elongated in 
some directions than in others, and in these cases the gravitational potential close to the 
horizon falls off more slowly than at asymptotically large distances. Nevertheless, if the 
limit I? — )• oo is taken keeping the horizon length scales finite, the region where this slower 
fall off occurs shrinks and disappears in the limit. We can confirm all these features by 
studying explicit known solutions. 

Myers-Perry black holes. Each independent rotation introduces a new length param- 
eter o. Even if this parameter is usually regarded as the ratio of angular momentum to 
mass, as we discussed earlier the length a is a more basic, purely geometrical quantity. We 
take these parameters a to scale in such a manner as to preserve finiteness of the metric. 

Let us take the limit -D — t- oo for Myers- Perry black holes with several non-zero rotation 
parameters. We find no important distinction between the even- or odd-dimensional cases, 
so for definiteness we take odd D. The solution is [26] 




{dt + ai^^dipif + 



n — r^r. 



liF 



dr^ , (2.19) 



where summation over i = 1, . . . 



D-l 
2 



is assumed, direction cosines satisfy fif = 1, and 



When ro = the metric 

ds'^ = -df + (r^ + aj){d^4 + /i- d^- ) + Fdr'^ (2.21) 

describes flat space in elHpsoidal coordinates. When ro > 0, the horizon is at r = th such 
that f'^rQ~^ = n(rj7). Let k be the number of non-zero rotation parameters Oj, and for 
simphcity of illustration set them all equal to a. Then we have 



1 + 72 ' (2-22) 



n \ r J V 
so rn is found by solving 

— =1 + ^ ■ 2.23 

If we keep k fixed when D — )• oo then rn — )• ro, and r^r|^^^/n clearly vanishes for r > rg: 
the metric becomes flat spacetime everywhere outside the horizon. 

If instead the number k of non-zero spins grows with D such that m = — k remains 
fixed, then in the limit — t- oo eq. (2.23) becomes 



So rn ro, but still 



2 \ 1/2 



+ ^ ^rH = Jrl-a\ (2.24) 



„2„n-3 /„2 , „2\ -0/2 



■ I 9 9 ; (2-25) 
n V + a 

vanishes, and thus leaves flat spacetime, at any r > rn when D ^ oo. This is the case 
even when ro ~ a, in which case rn is very small, the horizon can be highly pancaked 
(at least if some spins vanish), and rn *C a, which can be regarded as an ultraspinning 
limit [27]. The main point to note is that, in spite of the fact that at finite D there is a 
region very close to the black hole where the potential falls off like (instead of 

r^~^), the radial extent of this region vanishes in the limit D — t- oo when we keep ro and 
a fixed. 

Summarizing, in the limit — )• oo the gravitational field vanishes completely outside 
the horizon, and the effect of a 7^ is that the surface in (2.21) that is cut off at r = ro is 
not a sphere but an ellipsoid. 

The product GM of the solution vanishes in the limit for the same reason as in the 
static case. The angular momentum in gravitational units, GJ, also vanishes, which seems 
natural if we interpret that G vanishes. It is interesting that the ratio 

J 2a , , 

which is not suppressed by factors of J1d_2, approaches zero. On the other hand the 
angular velocity on the horizon 

= (2.27) 
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remains finite when D ^ oo, even though no rotational dragging is felt anywhere outside 
r = vh- The surface gravity in the singly-rotating case, 

' + 2.28 



2 \rjj + rn 
diverges like it did in the static case, but its minimum at 



a D-3 , , 

remains finite as -D — )■ oo. This suggests that the value a/rjj; = 1 separates the regimes in 
which rotating black holes at large D are stable or unstable to axisymmetric perturbations 
of the type discussed in [27] . 

Blackfolds. Other types of rotating black holes are known to exist in D > 5. Exact 
analytical solutions have only been found in five dimensions, but approximate solutions 
are known for black rings and many other black holes in arbitrary D > 5 [28, 29]. These 
are built by smoothly bending black branes, and are referred to as blackfolds. We can 
infer their large D limits from the properties of black p-branes. 

It is easy to see that a blackfold constructed out of a black p-brane, with p fixed 
as Z? — >• oo, becomes a 'dust-p-fold' as D — )■ oo when we keep finite the appropriate 
length parameters that determine its geometry {e.g., horizon thickness and worldvolume 
curvature radius). For instance, for a black ring with horizon x S^^^ we keep the 
radius of the S^~^ and of the fixed as — )• oo, and find a 'ring of dust' with vanishing 
rotation velocity and a fiat metric outside the singular horizon. 

We may also have black p-folds for which p diverges but n = D — p — 3 remains fixed 
as — )• oo. We are interested in the cases for which the gravitational potential beyond 
some distance Vc remains bounded above and behaves, at finite D, like {rc/r)^~^. Then 
when — 7- oo it will vanish for r > Tc. While it might seem possible that the radius Tc 
lies at a finite distance outside the horizon, the evidence we have found from the Myers- 
Perry solutions, which include blackfolds in this class, confirms our expectation that the 
spacetime becomes fiat outside the horizon. For these objects, GM must vanish. A simple 
example besides the Myers-Perry solutions are homogeneous black p-folds with volume 
VD-n-3, for which [29] 

IQttGM = {D- 2)VD-n-3^n+irJ^ . (2.30) 

We are keeping n and tq fixed. Now, in the case of e.g., a spherical blackfold, the volume 
R^~"'~^^D~n-3, and with it GM, vanishes due to Q^-n-s — ^ 0. So blackfolds of this type 
appear to conform to our general picture. There can also be blackfolds for which Vo-n-s 
does not vanish when L) — t- oo, for example toroidal blackfolds with Vd_„_3 = L^^^-^ 
with L fixed. Their mass-length (2.12) 

g-n-3 1 

^ r^-' {D - 2)T^ ^ L (2.31) 
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remains finite. It is doubtful tliat these objects can be sensibly defined in the limit D — )• oo: 
at large distances their metric coefficients cannot be finite. If a suitable limit existed, they 
would form a different, infinitely more massive sector of the theory, similar to the black 
p-branes discussed at the end of sec. 2.2. In this article we will not consider them any 
further. 

2.4 Charge 

The gravitational effect of the electric field in a Reissner-Nordstrom black hole is charac- 
terized by a charge-radius rg, related to the charge Q by 



This radius appears in the metric in the form (rQ/r)^(^^^\ If, as may be natural in 
this context, we assume that we keep rq fixed, then the gravitational effect of the charge 
vanishes at distances r > rq. Note that Q is measured in units of the gauge coupling (the 
electron charge). So if G is chosen to scale like (2.11) (to keep the mass fixed, for fixed 
ro), then we may keep both rg and the charge Q fixed if the gauge coupling is made to 
vanish like D~^/'^. 

Similarly, the effect of p-form charges (with p not scaling with D), which can be 
carried by black p-branes or by black holes with non-spherical topologies [30, 31, 32], 
vanishes in the directions transverse to the brane if we keep their charge-radius fixed. 
However, a p-form charge forbids the break up of the brane in the p directions parallel 
to the worldvolume. In these cases the equation of state of the brane does not become 
dust-like. For instance, for dilatonic p-branes in the limit D — )• oo one finds 

„ N sinh^ a 

P = ^e, 2.33 

l-FA^sinh^a 

where a is the charge-boost parameter (which remains fixed as D — )■ oo if the charge radius 
is fixed) and determines the dilaton coupling (see e.g., [32]). The field nevertheless 
vanishes outside the horizon. 

2.5 Large D in Anti-deSitter 

The cosmological constant A in Anti-deSitter spacetime introduces a curvature length 
scale. ^ In the limit — t- oo, instead of fixing A it is sensible to keep 



/ (0-l)(0-2) D 

^"V — w — m' ' ' 

since this is the length that appears in the metric. If we kept A fixed, we would have 
L — 7- oo and the effect of the cosmological constant would disappear from the geometry. 

''The first part of this discussion appUes also to de Sitter spacetime and to Schwarzschild-de Sitter black 
holes away from the Nariai limit. 
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Again, we focus on the simplest black hole solution, which has the same form as (2.1) 
but with [33] 



r 



/(r) = l-(^) +-^. (2.35) 



The event horizon is not at r = rg but rather at rjj > tq, the real positive radius where 
/(rn) = 0. At large D, 

= ro (l - ^ In (l + ^) + O {D-^)^ . (2.36) 

Thus, rn fo when D — t- oo. Like in flat space, in this limit the gravitational effect of the 
black hole outside this horizon vanishes. Black holes are again non-interacting particles, 
now in Anti-deSitter spacetime. However, although their area vanishes due to — ^ 0, 
these black holes are not always 'small AdS black holes', since rn can still be either smaller 
or larger than L. Their temperature is 

^_ {D- l)rj, + {D- 3)L^ 

so 



(2.37) 



Although this diverges when D — )• oo (like in flat spacetime), as a function of rg this 
temperature reaches a minimum at ro = L. We may then regard black holes as small or 
large according to whether ro < L or tq > L. We see in (2.36), (2.38) that the presence of 
L modifies the large- short scales, in particular the surface gravity length is 

and its dependence on rg is rather different in the two regimes ^ L. 

The shrinking effect of the Euclidean time circle and of the area of the spheres 
as — )• oo may be, if convenient, factored out by, e.g., dividing T by D, and dividing 
extensive thermodynamic quantities by For instance, the Euclidean action of the 

black hole spacetime vanishes, both because of the factor and because of the shrink- 

ing Euclidean time circle. In this case there is no trace of the Hawking-Page transition at 
D — 7- oo. However, if we introduce a rescaled free energy 



then the Hawking-Page transition at which F changes sign, is apparent at rn = L, al- 
though extremely abrupt: at D — )• oo this free energy changes from — oo for rn > L, 
to 0^ for rn < L. Another way to avoid these divergences is to consider ratios between 
extensive quantities. 

The limit of very large AdS black holes yields planar black branes in Poincare AdS. 
These black branes have a conformal equation of state P = e/{D — 1), so in the limit 
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D — )• oo they behave again like 'dust' branes, which are marginally stable to longitudinal 
fluctuations, with empty AdS spacetime outside them. In section 7.5 we briefly revisit 
these branes. 

In summary, all our knowledge about higher-dimensional black holes points to the 
conclusion that even when the efl'ects of rotation, different horizon topologies, charge, or 
cosmological constant, are accounted for, the limit D — )• oo consists of conflgurations of 
objects with flat (or AdS) spacetime outside their singular horizons. We return now to 
extract more information from the basic solution (2.1). 



3 Sphere of influence 

While in the limit D — )• cxd the gravitational fleld vanishes at all r > ro, when D is 
large but finite there is a small region around the horizon where the black hole exerts its 
gravitational infiuence. More precisely, outside the horizon the gravitational potential is 
still appreciable, i.e., 

(^)°-'=0(I,»). (3.1) 

within the region 

r-ro<g + 0(I)-2), (3.2) 

while it vanishes exponentially fast in D outside this range. 

For some observables this 'sphere of influence', characterizing the degree to which a 
phenomenon is localized close to the horizon, may extend out to a further radius. Typically 
this is due to radial derivatives of /(r), each of which brings in a factor ~ D. In general, 
the relevant result is that, for r > rg. 



D'i^-^) =0{D») ^ r-ro<j^{a + blnD) + OiD-'), (3.3) 

where a, b are D-independent numbers. For instance, the effects of the gravitational force 
drf have a range of this form with b = 1. We see this in the acceleration of an observer 
at constant r in (2.1), 

a=^^r-^Y-'^^r-^Y, (3.4) 

2rox/7(rO 2ro V r ; ' ^ ' 

which has a range 

r-ro<'^lnD + OiD-^). (3.5) 
An example with 6 = 2 is the Kretschmann scalar (2.7), which is of order one in the region 

r-ro<^lnD + 0{D-^). (3.6) 
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The scattering of light rays thrown towards the black hole (2.1) provides another 
illustration. These rays will be absorbed if the impact parameter b is smaller than the 
critical value [34] 

ro . (3.7) 



2 J \ D-3 

As — )• cxD this gives 6c — )• tq. Light rays that pass with b > tq suffer no deflection: it is 
infinitely difficult to catch a line of force that pulls towards the black hole. At finite but 
large D the photon sphere of influence is 

,,..„(l^l±M+0(I,-)). (3.8) 



The photon absorption cross section 



(3.9) 



vanishes as D^^/"^, as already discussed. For the rotating black hole (2.19), the critical 
impact parameter for light rays in the plane //^ = 1 asymptotes to be yj'^^o + ^fc ' which 
is precisely the radius of the circle at r = tq in this plane in the limiting flat geometry 
(2.21). 

3.1 Quasinormal modes 

Quasinormal modes characterize the black hole's own dynamics. It is natural to expect 
that they are localized close to the sphere of influence and controlled by the scale tq/D. 
This is actually the case. At large D the quasinormal frequencies can be estimated or 
computed (see [35, 36, 13] and section 6 below) with the result that for either scalar, 
vector or tensor perturbatons, 

D 

RecjQN ~ — , (3.10) 

at least when the angular momentum and overtone numbers are ^ D^. In the WKB 
approximation these quasinormal modes can be regarded as localized close to the circular 
photon orbit at 

1/(^-3) 
Tph = ( ) ^0 

The imaginary part of wqn is suppressed at large D relative to the real part, which 
may be related to the apparent ease with which black holes can fragment in this limit. 
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3.2 Across the horizon, a small interior 

In the limit — t- oo a particle that falls in the black hole encounters a singularity at 
r = ro where its trajectory comes to an end. However, for large but finite D the horizon 
is a non-singular place and we can ask how long it takes for the particle to reach the 
singularity at r = — i.e., what scale controls the internal size of the black hole. 

Taking for simplicity a particle that starts at rest at infinity in a radial trajectory, the 
proper time elapsed between the moment it crosses the radius r = R until it reaches the 
singularity is 

The two factors in this result each have a different origin. In the large D limit the particle 
takes a time that diverges exponentially with D to go from any finite distance outside the 
sphere of influence until it reaches this region, where (R/ro)^ ~ 1; this is just a reflection 
of the very small gravitational force outside this sphere. From the moment when the 
particle enters this region until it reaches the singularity, the time (3.12) that passes is 

t,,„3 ~ ^ + 0(D-2) , (3.13) 

i.e., a very short time, determined once again by the scale vq/D. Most of this time is 
spent in the sphere of influence, which now includes an inner region — r < r^jD. After 
that, the singularity is reached exponentially fast in D. Eq. (3.13) gives a sense in which 
we can regard the black hole interior as small. 



4 Classical radiation from black hole interactions 

The absence of a gravitational field inbetween two or more black holes implies that, when 
they move in the presence of one another, no gravitational radiation is emitted. Below we 
verify that this is consistent with the large D limit of several previously studied processes, 
but there is an effect that remains in phenomena that probe very short lengths < tq/D. 

4.1 Black hole collisions 

For simplicity we study collisions among black holes of equal size. It is clear that when D — >• 
CO two black holes thrown towards one another will not be deflected from their straight 
paths unless the impact parameter is < tq. An often studied case is the ultrarelativistic 
limit, i.e., the collision of two Aichelburg-Sexl shockwave solutions. It is immediate to see 
that at D — )• oo these solutions become trivially flat outside a singularity at the center of 
the shockwave plane — away from this point, the shock itself disappears. So, in this limit, 
for any collision course that is not exactly head-on, we conclude, unshockingly, that the 
particles simply fly by each other, with no emission of radiation. 
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Put now the black holes in a head-on collision course. Collision is unavoidable, but 
will there be any radiation emission? If the black holes start towards each other with 
infinitesimal velocities, then our argument about the entropy in a merger in sec. 2.2 tells 
us that no radiation will be produced when they merge. More generally, if they move 
towards each other with initial velocities if, then the area theorem imposes an upper 
bound on the ratio of the radiated energy to the initial energy [37] which, when D ^ oo, 
is^ 

e{v) = ^<l-y^l^. (4.1) 

As u — >• we recover e — t- 0, i.e., no radiation, but for ?; > there is the possibility that 
radiation is produced. This indicates that kinetic energy, but not any rest mass of the 
black holes, can be converted into radiation. 

We may expect that the actual e is maximized in ultrarelativistic collisions and in fact 
when — )• 1 the condition (4.1) imposes no constraint. For Shockwave collisions a more 
stringent bound on e can be obtained from the area of the apparent horizon at the moment 
that the Shockwaves meet [38, 12]. When D ^ oo this bound becomes 

Cshock < ^ • (4.2) 

According to refs. [18, 19], this bound is actually saturated, and thus gravitational radi- 
ation is indeed produced at L> — t- oo. This might appear to run against the picture of 
D — 7- oo black holes as non-interacting dust particles, but there is a plausible interpreta- 
tion for this result. A head-on collision is an extremely fine-tuned process (even more so 
at D — 7- oo where the cross section vanishes) which can excite dynamics at infinitely short 
distances involving modes of infinitely high frequency, such as the quasinormal modes with 
u ~ D/rQ, which are localized in the region (3.11) very near the black holes. At large D 
these modes can get excited if (and only if) very short scales ^ tq/D are probed. 

More precisely, the emission of radiation comes with a factor ~ oj^ from the frequency- 
volume available to radiation. If the characteristic length of modes that are radiated is 
i, then the process will be governed by a factor {uj£)^ . Emission will be completely shut 
off when the frequencies that can be excited in a process are to < but it may remain 
appreciable as D — t- oo when w ~ even if £ vanishes in this limit. In the following we 
find evidence that this effect is indeed present, and moreover that the relevant scale is 
i ~ tq/D, i.e., that of typical quasinormal modes. ^ 

Ref. [39] estimated the amount of radiation produced in black hole collisions in the 
'instantaneous collision approximation' (which may in fact become more accurate as D 
grows). At moderate (not ultrarelativistic) initial velocities, and ignoring D-independent 

^Even if the individual masses and areas vanish due to the shrinking factors of flD-2, these cancel out 
in this ratio. Also all factors of G cancel, so this is a quantity that can be defined with reference to purely 
geometric quantities. 

^The fact that the number of graviton polarizations grows like does not seem to play more than a 
subleading a role, negligible for the previous argument. 
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numbers, the leading large D emission per solid angle is^ 

^ ~ GM^u:^' , (4.3) 

where is a (physical) cutoff in the frequencies that are radiated. Integrating over all 
angular directions, the ratio of the total radiated energy to the black holes' mass is 

^^1m m dn' ^ \ D ) 

(we have used (2.12)). If the frequency cutoff were ujm ~ then no radiation would 
be emitted when D — )• oo. But precisely at the frequencies Um ~ D/vq there is a big 
enhancement in the emission of radiation. 

Our next example provides another instance in which there is a possibility of producing 
radiation by exciting frequencies ^ D/tq. It is worth stressing that neither of these 
calculations involve a quasinormal mode study, nor indeed any perturbation analysis of 
the black hole spacetime (2.1). We find remarkable the consistent appearance of this scale. 

4.2 Orbiting black holes 

The radiating power from two black holes, each with Schwarzschild radius ro and mass 
M, orbiting around each other at a distance / and with orbital frequency cj, tends at large 
D to s 

~ Qd-2GMH^u:''+^ (4.5) 

(again up to D-independent numerical factors) [39]. The fraction of the system's energy 
emitted per orbit is 

If the black holes follow Keplerian orbits with I > tq then 



^^-1 I \Vdi 

vanishes as D grows, since their mutual attraction becomes increasingly weaker (this 
also justifies the simplifications made in deriving (4.6), although note that the orbits are 
unstable). This does fit the picture of non-interacting, non-radiating particles. 

However, let us push the argument further and consider that the black holes follow 
circular trajectories under an external force (whose gravitational effect we neglect; we 
proceed even it is unclear how sensible this is), so that u is independent of the black hole 
masses and their separation. The dominant factor in (4.6) at large D is {utq/ D)^ . So 
when the orbital frequency is a; ~ E)/rQ, i.e., when the motion is so fast that it is capable 



'^The dependence on dimensionful parameters in this equation can be easily worked out from generic 
considerations. But note that we are also accounting for the dominant D-dependent dimensionless factors. 
*See footnote 7. 
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of exciting the quasinormal modes of the black hole, radiation is possible as D — )• oo. In 
this regime the approximations made in deriving (4.6) certainly cease to be justified, but 
still we can observe how the excitation of the very high-frequency modes in the thin region 
near the horizon may result in radiation that survives the limit — )• oo. 

The themes of this section will reappear in sec. 6, where we confirm this picture of the 
interaction of the black hole with classical waves. 



5 Quantum effects 

We do not study quantum gravitational effects in any detail but merely make some ele- 
mentary remarks. 

Since the curvature (2.7) at the horizon grows unbounded at large D, one can expect 
a breakdown of the semiclassical gravitational effective theory in the close vicinity of the 
horizon. The same effect is seen in the black hole temperature, 

As D — 7- oo, Th reaches arbitrarily high values eventually hitting the Planck scale. This 
occurs even when tq is larger than the Planck length Lpianck; because the scale controlling 
the curvature and temperature is not ro but vq/ D. In this reasoning we are assuming that 
the Planck scale does not scale with D. If, instead, we choose to scale it as -Lpianck ~ 
(and so -Epianck ~ D), e.g., Lpianck = -^Pianck/-D, then quantum gravitational effects 
might be suppressed when tq » -C-pianck- This type of scaling, therefore, might seem more 
appropriate for the study of semiclassical effects. 

The situation, however, is subtle. Consider Hawking radiation at large Z?, which was 
studied in ref. [16]. One might expect that the typical frequency of Hawking quanta 
should be of the order of Th ~ D /r^. However, the actual typical frequencies are much 
larger. The reason is the huge increase in the phase space available to high-frequency 
quanta at large D (already encountered in sec. 4), which shifts the radiation spectrum 
towards frequencies much larger than Th- The number density distribution of quanta at 
temperature T is 

n{uj)duj = VLD-2—-pf. — T, (5.2) 

with the factor coming from the angular integration in momentum-space. At fixed 

but large this density peaks at [16] 

a; = L>T(l + 0(e~^)) >T. (5.3) 

Then, for a black hole the typical frequency of Hawking quanta is 

UJH ^ DTh ^ . (5.4) 
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We see that even if we scaled -Epianck ~ -D, we would still have w// /E'pianck ~ D ^ 1: 
Hawking quanta would still reach superPlanckian energies. It is unclear whether a stronger 
scaling of the Planck length with D would be useful. 

Thus the black hole is a very large quantum radiator, whose radius (~ tq) and typical 
classical vibrational wavelengths (~ vq/D) are much longer than the wavelengths it radi- 
ates quantum mechanically (~ tq/D'^). In this respect it is interesting that, even if phase 
space integrals receive suppressing factors of from angular integrals (from both mo- 

mentum and position space), these are offset by the growth ~ uj^~'^ of volume in radial 
(frequency) directions. The number density (5.2) at the peak frequency (5.4) behaves like 

Then, the total Hawking flux, even if radiated by an object of vanishingly small area 
Ah ~ D^^/^, still shows factorial growth ~ . This enhancement is due to the fact 
that Th ~ Z).^ For practical applications of the large D limit, the ultrashort wavelength of 
Hawking radiation can be a bonus since it implies that the geometric optics approximation 
for graybody factors applies very accurately [16]. 

It would be interesting to investigate whether the black hole information problem can 
be sensibly formulated in the l/Z? expansion. Let us mention, among the potentially 
relevant factors, that at — t- oo the singular horizon (which has zero area) perfectly 
reflects radiation of any finite frequency. This might seem to prevent information from 
being lost across the horizon (although it might also be destroyed there), but on the other 
hand Hawking pairs come out at infinite frequency. For very large but finite D the very 
short time (3.12) available for measurement for observers who cross the horizon might 
facilitate the consistency of a form of black hole complementarity, but the extremely fast 
evaporation rate compresses enormously all the timescales involved. It is then unclear 
without a more careful investigation whether these are positive or negative features for 
making this a useful approach. 

6 Large D effective theory: scalar wave absorption 

Up to this point we have been mostly drawing consequences from the limit D — )• oo in 
processes that had already been calculated at finite D. We have gathered ample evidence 
that a consistent picture emerges. Now we apply the large D expansion to problems that 
have not been fully solved previously in analytic form. The main idea that simplifies their 
study is the following. 

6.1 Near and far regions 

The gravitational field is strongly localized within a region close to the horizon character- 
ized by the scale vq/ D. This is a length that at large D is widely separate from the radius 
®See footnote 6. 
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ro, so we can define two distinct but overlapping regions in tlie geometry: 



near region : r — tq <C tq , 

far region : r — rg S> — . (6-1) 
Equivalently, introducing the variable R = (r/ro)^ they can be characterized as 

near region : InR <^ D , 

far region : In R ^> 1 . (6-2) 

The far region, where the geometry is effectively flat, excludes the sphere of influence (3.2). 
The latter is part of the near region, ^'^ which extends into an 

overlap region : — <^ r — <^ vq , i.e., l^lnR^D, (6-3) 

very close to the horizon, where it smoothly connects to the far region. Then, the study 
of any phenomenon that takes place in this geometry lends itself naturally to the method 
of matched asymptotic expansions (first used in this context in [7] ) . If we manage to solve 
the field equations in the near region, with regularity conditions imposed on the horizon, 
then this solution will provide boundary conditions for the far field by requiring that the 
fields match where they overlap. 

6.2 Massless scalar wave equation 

We study massless scalar field propagation in the black hole background (2.1). Throughout 
this section it will be slightly convenient to work with the parameter 

n = D-3 (6.4) 

instead of D. At large D they are of course equivalent. 
We study the wave equation 

D^- = (6.5) 

in the background (2.1). We set 

^ = e— *V./(r-)5;!i(0), (6.6) 

where l^^^^(r2) are spherical harmonics on S^'^^, and henceforth we drop the mode indices 
from ip. The equation becomes 



1 d 



It is conventional to write this as an equation for 

(j){r) = r("+^)/2^X^) 



"incidentally, the near region retains more features than a Rindler horizon. 
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in the form 

+ = 0, (6.9) 

where the tortoise coordinate defined by dr^, = dr/f{r) is given in terms of a hypergeo- 
metric function, 

r.=.F,f-i,l,^;n"V, (6.10) 



n n \ r ^ 
and 

nn)={^((2/ + n)^-l + (n+l)^(^)") (6.11) 

is the Regge- Wheeler potential for these perturbations. 

At large n this potential scales like n^. In order to capture physics of interest we 
introduce 

Lo=-, l = - 6.12 
n n 

and consider uj and / as 0(1) quantities. Now when n is large the equation becomes 

^^+n'[6j'-V{r,))cP = 0, (6.13) 

with 

V(r.)='^{(2Ul?^{^y). (6.14) 

This potential vanishes on the horizon at r* — )> — oo and in the asymptotic region at 
— >• +00. It has a maximum at 



ro(2n^|±il + 0(n°)l . (6.15) 



When I » 1 this is 

W = ro (l + ^^^) + 0(n-2) , (6.16) 

which reproduces the radius of the circular photon orbit (3.11). The maximum of the 
potential, corresponding to critical scattering at the threshold of absorption, occurs for 
the frequency to = Uc, with 

uJcro = i+^ + 0{n-^). (6.17) 

This gives the real part of the quasinormal mode frequency (3.10) in the WKB approxi- 
mation. The critical impact parameter be = l/ojc at / ^ 1 reproduces the geometric optics 
result (3.9) at large n. This critical frequency will play a central role in the problem. 

When n — y oo this potential becomes extremely simple: for r^< > ?'max (where n = r) 
we recover the flat space potential with only a centrifugal barrier. For r^, < rmax the 
potential vanishes exponentially quickly in n. So the limiting form of the potential is 

V{n)^'^e{n-ro). (6.18) 
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We recognize here two of the main running themes of this article. First, waves outside 
the black hole propagate in a flat space potential. Its height F(rmax) ~ "-^ becomes infinite 
when n — )• c«, so in this limit excitations of frequency co = 0(jiP) are not absorbed at all: 
this radiation does not interact with the black hole. Second, waves of frequency to > n/rQ 
can penetrate the barrier and probe the dynamics near the horizon of the black hole. 

Problems of field propagation are characterized by the field amplitudes at infinity and 
at the horizon, 

' (e--'-* + Ri{u)e'^-') r, ^ oo , 

Af{u)Ti{u)e-^^''* n^-oo. 

When the amplitude of the incoming wave Af^{uj) is nonzero, the physical information 
is contained in the reflection and transmission {i.e., absorption) amplitudes, Ri{u) and 
Ti{uj). We will solve this problem of scattering and absorption by the black hole, and in 
particular compute the absorption probability 

ji{u;) = \Tiiu;)\\ (6.20) 

The condition A™{u) = defines a different kind of problem: the determination of the 
spectrum of quasinormal modes of the black hole, which appear as poles in Ti{uj) and 
Ri{oj). This is a subtle eigenvalue problem to solve in the near region, and we postpone 
its detailed study to future work. 

6.3 Integrating out the near region 

In order to simplify the notation, we now set 

ro = 1, (6.21) 

so that 

/(r) = 1 - r"" . (6.22) 

The tortoise coordinate r^: is not quite appropriate in the near region — oo < r^: < rmax: 
the non-trivial features of the potential (6.14) in this region are erased in the limit (6.18) 
in which r^: remains fixed as n — )• oo. Although we could instead keep r^, = nr^, fixed, we 
find more useful to work with the coordinate 

R = r'^, (6.23) 

in terms of which the wave equation (6.7) becomes 



To leading order at large n in the near region, where In R ^ n, this equation becomes 
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where instead of / we use Wc, defined in (6.17). Tlie expansion breaks down when uj is of 
order n or higher. Then our calculations apply in the range w <C n^. 
The general solution to this equation is 



V^(R) = ^i(R - 1)-*'^ 2Fi(g+, + 1 - R) 

+yl2(R - 1)*" 2^1(1 -q+A-q-,2-q+-q.-l- R), 



where 



The regularity (ingoing) boundary condition at the horizon (6.19) takes the form 



V'(R) 



R=l 



OC 



(R-1)-^^ (l + 0(n-i)) 



(6.26) 
(6.27) 

(6.28) 



and therefore requires that ^42 = 0. We use the arbitrariness in the overall amplitude to 
set Ai = \l\fn for later convenience. 
The solution can now be written as 



V.(R) = ^('^+ + ^-^ (R-l)- 



n 



r(g+-9-) 



+R-«+ 



r(g+)2 
r((?- - g+) 



2Fi(g_,g_,l-g+ + g_;l/R) 



2i^i(g+,g+,l + g+ -g-;l/R) 



r(g_)2 

Expanding at large R yields the solution in the overlap region 1 ^ In R <C 



(6.29) 



V^(R) 



r(g+ + g_) /r(g+ - g_) 



-R 



n 



r(94 



-R'?+ + 



r(g-- g+), 
r('z-) 



R^- ) +0(R-^). (6.30) 



When (7+ = g_, i.e., a) = Wc, this expansion is not valid, and instead one gets terms 
R-i/2 R. We expect that this is due to the presence of quasinormal modes. Logarithmic 
terms also appear when — g_ € N. 

The flux of the scalar at the horizon, derived from the wave equation (6.9), is 



■A. 



m , 



)— 
d 



d 



dR 



dR 



R=l 



For our solution (6.29) we find 



horizon 



(6.31) 



(6.32) 



Eq. (6.30) is an important result in this analysis: by providing a boundary condition 
for the fields that propagate outside the 'sphere of influence', it codifies the physics of the 
region where all the black hole dynamics is concentrated. 
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6.4 Far region waves 



In the far region we set /(r) — t- 1 + 0{e the wave propagates effectively in flat space, 
obeying the equation 



1 d 



This is solved in terms of Bessel functions, 



V'(r) = 0, 



(6.33) 



r 



n/2 



rn/2 



From the behavior at large r 



1 



(Ci - iC2)e''^' + (Ci + iC2)e- 



V27ra; r"+i 

we infer the incoming amplitude and the incoming flux from infinity. 



1 

2tt 



Ci + iC2 



(6.34) 



(6.35) 



(6.36) 



Quasinormal modes would be obtained under the condition Ci + iC2 = 0. 

In the overlap region, and to leading order at large n, this solution gives (see ap- 
pendix B.l) 



V^(r) 



l + q++q- 

R ^ 



V 2TTnoJc tanh uq 



with 



, \j2T^nLiJc tan /3o 



1C 



{Ci - iC2)K^ R'^ + ^'t'^^ R'- 

ri-,7, 



CO < Uc, 



CO > LOc, 



(6.37) 



-nuJc (cfo —tanh ao ) 



CO < LOc , 



-mcJc(/3o-tan/3o)-i7r/4 ti) > 



(6.38) 



and ao and /3o defined by 



CO 
COr. 



sech ao , < Wc , 
sec /3o , < w , 



so that 



q+ — q- = 2\J co"^ — Qo^ = 2lOc tanh ao = 2icOc tan (3q 



(6.39) 



(6.40) 
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6.5 Matching and analysis of results 

Matching the coefficients in eqs. (6.30) and (6.37) we find 

^/■K{q+ - q-) T{q+ + q_) T{q+ - 



Ci 



Kc,T{q^ 



Kco^jT^{q+ - q -)T{q+ + q-)T{q- - q+) 
2T{q. 



C2 = - -^-'^^^ '"ow J ^ (6-41) 



when (2} < ojc, and 



Ci - iC2 



Kc,T{q, 



n i n Kcy/m{q+ - q^)T{q+ + q-)T{q- - q+) 

Ci+iL2 = T{q_Y ■ ^ ' 

when oj > uJc- 

This solves the problem of scalar field propagation in the presence of the black hole, 
since using these results we obtain the reflection and transmission amplitudes 

Ci-iC2 , . . . 

which satisfy + jT/p = 1. A simple quantity of interest is the absorption probability, 

/ V -^horizon ll^Gj i n a a\ 

-n{uj) = — = — (6.44) 



^in \Ci+iC2\ 

Observe that all the dependence on n in Ci and C2 is contained in the factor K^j. 
6.5.1 Low frequency: Co < Uc 

Since oq — tanhoo > 0, in this regime -R'^ is exponentially small in n, except for Uc — uj 
where Ki^j rapidly approaches 1. We can write 



^!l±^MEE.\ . (6.45) 

Given that Kijj <^ 1, the amplitude is dominated by \Ci\ ^ IC2I and we approximate 

2-11(2; 



„ (6.46) 



Vw2_^2r(^^ _ g_)2 |r(i + i2u)\ 

which is strongly suppressed. More explicitly, if ^ cjc, and since ojc > 1/2, we can write 
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with all the other factors in (6.46) remaining of order one. Thus, restoring the radius ro 
and oj = noj for clearer illustration, we conclude that waves of frequencies oj < n/rg are 
very strongly reflected by the black hole and interact very little with it. 

At very low frequencies oovq <C 1 we can check against earlier results. Consider s- waves, 
I = 0, i.e., ujc = 1/2, which are the dominant component of the absorption. In this case 
we have 

7K-) ^ ^-^Kl (6.48) 



,2 /ewroV" ' "2-2 



with 

where we have used Stirling's formula. Introducing the horizon 
find 

7(-).^o = (l - ^ + 0(a;3, n^)) . (6.51) 
Finally projecting plane waves onto s-waves, we find the scalar absorption cross section 

<br' = Ah(i-"^ + 0(a;^n-2)) . (6.52) 
The leading term is the universal result of [40]. 

6.5.2 High frequency: Cj > 

In this case Kcjyu)c is purely imaginary, with 

\K^>u.^ = l. (6.53) 

When cD S> (^c it can be approximated by K^j^^^ ~ e^*""^. 

In contrast to the previous regime, the absorption probability now does not depend on 
Kcj. It takes the simple form 

sinh(27rci)) sinh ( 2'K\JCj^ — uj^ ) 

7K^) = — ^ ZTzi ■ (6.54) 



^cosh vr ^d) + \J Qp- — 



When cD S> (^c this becomes 

7/(u;) = 1-0 , (6.55) 

which is the expected result: the black hole is an almost perfect absorber at very high 
frequencies. 

Note that the absence of a dependence on n in (6.54) implies a scaling behavior with 
Lij and I of the absorption probability at large n. 
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6.6 Effective scalar wave dynamics 

Matching constructions like we have performed admit an interpretation as effective the- 
ories. In our case case, the theory consists of scalar waves that propagate in flat space. 
When they reach a horizon (at r = 1) the field modes tp^iir) must satisfy the boundary 
condition (from (6.30)) that 



LUl 



n(, r(g+ - q-)T{q. f - r(g_ - q+)T{q+f 

2 - '-^nq, - q^nq.y + r(,_ - ,,)r(,,)^ - ^ ) • ^'-''^ 



This equation encodes all the scalar dynamics of the black hole, to leading order at large 
n. Once this condition is imposed, the reflection and transmission amplitudes for field 
propagation take the form that we determined above. 

This effective description differs in important respects from the one that results from 
the more familiar study of black hole absorption at very low frequency oj <^ l/rg [41, 42, 
43]. That analysis also performs a matched asymptotic expansion between a near region 
r <C UJ~^ and a far region r ^> tq. Then the effective theory is obtained after integrating 
out the degrees of freedom at scales < tq. In the large n effective theory, instead, we 
integrate the physics near the horizon at scales < vq/D. 

This distinction is important. The conventional low- frequency effective theory treats 
the black hole as a point particle: waves with uj <^ 1/ro do not resolve the size of the 
horizon. The conditions on the far field are then effectively imposed at r = 0. In our 
large n effective theory, instead, the horizon radius ro = 1 remains finite (even if the 
horizon area vanishes at n — t- oo), so the boundary conditions are imposed on a sphere of 
finite radius. This is why we are able to capture a much larger frequency range, including 
CO > n/ro ^ l/?'0) in which the black hole acts as an almost perfect absorber: excitations 
of these wavelengths perceive it as a finite-size object. 

6.7 Final remarks 

Upside. We have obtained a solution for the scattering and absorption problem that 
gives, in a simple analytic expression, the amplitudes over a very wide range of frequencies 
(and partial waves), running from the lowest part of the spectrum until very high values. 
As an illustration, we not only recover the well-known universal result for the low-energy 
absorption cross section but also we can easily extract corrections to it, eq. (6.52), which 
would be very hard to obtain in other approaches. 

The results also confirm the conclusions of section 4: the interaction of the black hole 
with waves of frequencies that remain fixed as D grows is strongly suppressed by a factor 
of the form ~ (uro/D)^ (see e.g., eq. (6.47) or (6.49)), while frequencies that scale like 
D interact appreciably with the black hole. In fact waves with 

+ (6.57) 
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are almost perfectly absorbed. 

Our analysis also provides the solution for several types of linearized gravitational 
perturbations of this black hole [44]: perturbations that are tensors on 5*^"^ are governed 
by the same equation as (6.7) at all D, and scalar gravitational perturbations also obey, to 
leading order at large D, the equations (6.13), (6.14) when / = 0{D^). Moreover, gauge 
field perturbations satisfy the same equations at large D including I = 0{D^). Shear 
perturbations of black branes are also known to obey the equations for massless scalars. 

Downside. During the matching construction we have found that it breaks down in 
some specific instances. 

The first one occurs when — q~, i-e., when uj ~ ujc- This is also the region in which 
IC1/C2I — 1- Since the quasinormal modes appear when Ci + iC2 = 0, we see that our 
present construction is not well suited for the calculation of their frequencies. Nevertheless, 
this analysis shows that quasinormal modes must be expected when w ~ cjc. This agrees 
with the WKB considerations in sec. 6.2 and refs. [35, 36, 13]. Quasinormal modes are of 
course very important since they are responsible for resonant scattering which results in 
damped ringing at their frequencies. We expect that a modification of our technique will 
allow an analytic solution of the quasinormal modes and their frequencies. 

Second, our matching is also incorrectly performed when — q- G N. While the 
matching result for Ci may still be correct, the matching of C2 is more complicated. 

Our solution should lose accuracy around these special frequencies. For instance, at 
all of them the absorption appears to vanish since C2 diverges, but this conclusion is not 
reliable. In constrast, for all uj > ojc there is no problem in matching Ci and C2 and our 
results in this range away from ~ Wc are sound. 

Finally, the effective theory approach fails to apply for the ultra-high frequencies and 
angular momenta in the range ujvq, I > D^. The wavelength of these excitations is so short 
that they are insensitive to the curvature and they do not distinguish between near and far 
regions. Still, our expression for the absorption probability can be expected to smoothly 
merge with the result at these ultra- high frequencies when coro > D^, I, since in this case 
the deviations from 7;(ci;) = 1 must be extremely small. However, the range euro ~ I ^ 
seems to lie outside the large D techniques of this section. 

Aside: near-horizon conformal symmetry? The appearance of hypergeometric 
functions in the near-region solution is suggestive of the presence of a two-dimensional 
conformal symmetry governing the amplitudes [45, 46]. Although in this article we will not 
pursue this idea, let us mention a potentially relevant fact that arises when g+ — g_ G N,^^ 
which is always in the regime to < ujc- For odd q^ — q^ = 2m+l, the absorption probability 

^'^In these cases the matching for C2 is not vahd, but we expect that the resuh for Ci, which is what we 
use, remains valid. 
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(6.60) 



the absorption factors can be written as bosonic or fermionic Boltzmann thermal functions 
at temperature Th (or two sectors at Th/2). 

The conformal symmetry in [46] appears in the low frequency regime oor^ ^1. As 
discussed in sec. 6.6, our approach can deal with much larger frequencies, capturing the 
dynamics up to a; < D^/ro. Then the amplitudes and absorption probabilities above, 
although similar to the ones in [46], differ from them. It would be very interesting if the 
large D effective description for the black hole took the the form of a two-dimensional 
conformal theory. 



7 Black brane instability 

In the previous section we have solved a scattering problem, namely the interaction of the 
black hole with waves that propagate outside it. In this section we investigate a problem 
closer in spirit to quasinormal mode analysis, where we study the dynamics of the black 
hole itself. This is, we solve an eigenvalue problem for the field in the near region with 
boundary conditions that correspond to the absence of any external perturbing sources 
outside the black hole. The modes we seek differ from quasinormal ones in that they 
correspond to an instability. 

Refs. [23, 24] proved, through a numerical solution of linearized perturbation equations, 
that the black branes of (2.15) are unstable to fluctuations along their worldvolume, 
^9^lv ~ gf^*+«k z^ Besides demonstrating the large D expansion in a different problem, 
the determination of the spectrum ^[k) is an excellent benchmark for large D studies: 
good numerical results are available in several dimensions, but also there exist analytical 
approximate results [6, 7, 9, 21]. Thus we can test both the accuracy of the method when 
applied to finite values of D, and its effectiveness in comparison to other approaches. The 
outcome is very good on both counts. 
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7.1 Perturbation equations 

At linearized order the perturbation problem depends on 

n = D-p-'i (7.1) 

but not on D and p separately, so in this section we will use n as the expansion parameter. 
To lighten the notation we set tq = 1 like in the previous section. 

The Gregory-Laflamme problem involves perturbations that are scalars on S""*"^. There 
is only one physical degree of freedom, which by an appropriate choice of gauge can be 
chosen to be 

8gtr = r/(r)e™+*'=^ , (7.2) 

all other metric components being obtainable from this one. The linearized perturbation 
equation is 

r]" + P{r)r]' + Q{r)r] = (7.3) 
where the functions P(r) and Q{r) are [24] 

P(r) = ] \ 3n^ - UnQ'^r^ + {3n'^ - 6n^ - Snk'^r^ - Ar'^Q^ + Snfl'^r^)/ 



rfA{r) 



{6n^ - 3n^ + 4k^r^ - Ank^r'^)f + f , (7.4) 



+ (3n3 - + r?k^r'^ + lOn^f^V^ + 8k^n'^r^)f 

+ {n^{l - 6n - n^) + 2k^r'^{2k^r^ - 2n - n^) + n^r^{4: + An - Sn^)) f 
+ i-2n^ + 3n^ + + Ak'^r^ + SnfcV^ + n'^k'^r^)f + n^/^j , (7.5) 

with 

A{r) = 1? - 40^2 - (4/cV2 + 2-r?)f + f . (7.6) 

We shall solve the boundary value problem that results from requiring asymptotic 
flatness at infinity and regularity at the horizon. Note that P and Q have a simple pole 
at r = Ts > 1 where A{rs) = 0. This is a regular singular point of the equation. We have 
analyzed this point and checked that the solution that we find below is indeed regular 
there. 

Taking cue from the result (2.18) that the zero- mode wavenumber scales like l/^/n^ 
we define 

k = ^, (7.7) 



and regard k as being 0(n*^), which allows us to keep track of very short wavelengths 
~ r^^l^Jn. In contrast we take 17 = 0(v^\ 

We solve the problem by matching the solutions in the near and far regions of sec. 6.1. 
Our method is an adaptation and extension of ref. [7]. 
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7.2 Far region 

In this region r~" is a quantity that is exponentiaUy smah in 1 /n and therefore does not 
yield any perturbative 1/n corrections. Eq. (7.3) for r]'^^^'^^ becomes the flat space equation 



dr^ r dr \ r"^ 

Introducing the notations v = (n + 2)/2 and /c^ = \/ nk'^ + Q'^, the solution that is regular 
at infinity is 

where K^{x) is the modified Bessel function of the second kind. We fix the integration 
constant for the amplitude of the perturbation in a manner that will simplify the matching 
in the overlap region. Here we use the coordinate R defined in (6.23) and expand up to 
next-to-next-to-leading order, 

InR (lnR)2 ^, , , 

r = 1 + + ^—f- + 0{n~^) . 7.10 

Expanding now (7.9) for large n we find 

-MW ^ 4 / 1 1 + P In R 
W2 R n R 



^J^2(.^-l.-)lnR+(l + .^)(lnR)^^ j_ (7.11) 

2R ' 



Here is a factor (independent of R) in which we absorb all the n-dependence from terms 
of the form oc l/(n*R) in this expansion. We choose the integration constant in rj^^^^^ to 
eliminate this factor, so that 

ffar) _1 1 + fc^lnR 

' R n R 

12(^^-n'^)l„R + (l+^^)(lnRy^ + 0(„-). (7^12) 
2R 

Two comments about this result. First, note the absence of a constant term oc R*^. This 
is a consequence of asymptotic flatness. The second independent solution to (7.8) is 
Iuiknf)/f'"'^'^ , which grows at r — t- oo: this is a non-normalizable perturbation, which 
amounts to introducing sources for the field at infinity. This second solution would also 
yield a term oc R^ in the overlap region. Therefore the requirement that in this region 

lim r] = 0, (7.13) 

R— >oo 

is equivalent to asymptotic flatness in the far region. 

Second, the expansion (7.12) involves higher powers of InR but not of 1/R (this is 
true also if we include non-normalizable perturbations), essentially owing to the strong 
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localization that makes / = 1 + 0(e~"). Nevertheless, it is possible to compute the 
terms that come from the next-order far solution. Although we will proceed without them, 
they may be used to provide an alternative matching calculation. In appendix B.2 we give 
some details of this. 

7.3 Near region 

We expand in 1/n as 

"-E^^ (7.4) 

Eq. (7.3) gives an equation at each perturbative order of the form 



d 



R(R - 1)'^) + (2R - i)(R - 1)%) = ^o), (7.15) 



where is a j-th order source term built out of the solution up to order j — 1. The 
homogeneous equation has as its independent solutions 

1 ln(R-l)-lnR , , 

With these the solution to eq. (7.15) can be found using Green's method in the form 
r?Q) = AjUo + BjVQ 

+no(R) voiR')S^,)iR')dR' + voiR) no(R')5(,)(R')dR'. (7.17) 

The integration constants Aj and Bj will be determined by matching to the far region 
solution (7.9) at large R, and requiring regularity on the horizon at R = 1. This last 
condition can be derived by solving eq. (7.3) near r = 1. The regular solution is 

rj oc (r- l)-i+^/"(l + 0(r- 1)) 

+^(ln (R - l) f + 0(n-^ R - 1) j . (7.18) 

Matching condition. The overlap region corresponds in the near region to R ^ 1. 
Eq. (7.13) then implies that 



limy^ = 0. (7.19) 

j 



This condition can easily be seen to imply that the Wronskians of uq and i^q-j , 

W[uo,VU)] = R(R - l)'N(R)r?'y)(R) - n'o(R)r?(,)(R)] , (7.20) 
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satisfy 

^^^1 WWmI^ (7.21) 

j 

Since these Wronskians are given by 

W[uo,V(j)] = j dRuoS^j), (7.22) 



then the boundary condition in the overlap region can be conveniently written in the form 

R->c 



,l-^E;^/'^R-o5(,)=0. (7.23) 



Crucially, this j-th order boundary condition can be imposed with knowledge of the solu- 
tion just up to (j — l)-th order. 

We can now proceed to solve (7.15) order by order, with the sources 5q) computed in 
appendix B.3 and with the boundary conditions (7.18) and (7.23). 

Zeroth order. The equation at this order is homogeneous, so 

Ao ln(R-l)-lnR , , 

Horizon regularity (7.18) fixes Bq = 0. Matching to the far solution in the overlap region 
(7.12) fixes the amplitude Aq = \. 

First order. The condition (7.23) for r/(i) is automatically satisfied. We can integrate 
the source terms in (7.17) and determine the constants A\, B\ by imposing the boundary 
condition at the horizon and matching with the far region solution. Thus we obtain the 
first order solution in the near region as 

l)ln(R-l)-(l + P + r^)lnR 
^(1) = r31 • (^-25) 

Second order. The condition (7.23) to this order requires that 

-^(rj2-P + 2Prj + A;^) = 0. (7.26) 



Choosing A; > 0, the solution to this equation that gives unstable modes is 

n = k-P. (7.27) 

This gives the dispersion relation for the Gregory-Laflamme instability to leading order at 
large n. 
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As discussed above, this result does not require the actual second-order solution: only 
the source term computed with rj^i-^ has been used. However, we will need 77(2) ™ order to 
find further corrections. Integrating eq. (7.15) for j = 2 we find 



V{2) 



6(R-1) 

+ ^ — (R - l)f + (20 + 2P + l)(ln R)^ 



2fR-ll 



+2(fc^ - 17^ - n{l + P + n) In (R - 1)) In R 
-Aik^ + + 2Pn-P){R-l) 

-2{n^ + k^ + 2Pn - 2p)Li2(l - R)) , (7.28) 



where Li2(^;) is the dilogarithm function. The integration constants A2 and B2 have been 
fixed by matching to the far solution (7.12) and to the horizon solution (7.18). In more 
detail, near the horizon the solution becomes 



l) + ^(ln(R-l))2 + 0(R-l)], (7.29) 









R=l 












n 



where 



B = . (7.30) 



Third order. The regularity condition (7.23) for 7/(3) becomes 



2 



n '^ L 



f]2 _^ 21^3 _^ 2p - 2Pn + 2PJ72 + A:^ - 2k'^n - 2k^ 
+ (1 + P)(172 - P + 2Pl7 + P) In R 
= 0(n-4), (7.31) 

whose solution 

k ^ , 

Q. = k-k^ (1 + 2k- 2P). (7.32) 

2n 

gives the dispersion relation to next-to-leading order in 1/n. Again, this has been obtained 
before calculating 77(3) from (7.17). 

The boundary condition at the horizon is satisfied by setting ^3 = BVl in (7.30). The 
second integral in (7.17) can be performed analytically. This implies that we can compute 
the behavior of the third order solution 77(3) in the overlap region explicitly, which is 
crucial in order to impose the regularity condition (see appendix B.3). The constant A-^ 
is determined by matching to the far region solution (7.12). Although this can be done 
explicitly, it turns out that ^3 does not enter in the condition (7.23) for the fourth order 
solution. So we can proceed to the next step without specifying A^. 
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Fourth order. With the previous solution we can compute the source 5(4). Condition 
(7.23) at fourth order gives 

2n ' 

k . . ^ , . 

+—^(9 + 2Ak + 12/c2 - SiT^k^ + Svr^A;^ _ i2A;4). (7.33) 
24n^ 

In order to proceed beyond this point we should do the first integral in eq. (7.17) to 
impose the boundary condition on 77(4) at the horizon. However, we have not managed to 
do this analytically. Then we cannot determine -B4, which affects the fifth order regularity 
condition that would give the 1/n^ term in Q{k). So we stop at this order. 

Eq. (7.33) is the main result of this section. 

7.4 Comparisons and accuracy 

Setting J7 = in (7.33) gives the wavenumber k of the threshold zero-mode. Reverting to 
k = ^/nk^ we find 

^CL = ^/^(l-^ + ^ + 0(n-3)). (7.34) 

This reproduces the result in [7], which was obtained with a method essentially similar 
to ours, but using a different gauge. We have discussed the interpretation of the leading 
order result earlier in (2.18). 

Analytic approximations to r2(fc) have been computed from a rather different approach. 
Refs. [9, 21] solved black brane perturbations in a hydrodynamic expansion at small k for 
arbitrary n. Ref. [9] conjectured that the relation J7 = A; — /c^ is exact when n — )• 00. Our 
results, already in (7.27), do prove it. Ref. [21] extended the calculation to include terms 
up to oc k?' . These hydrodynamic results and our large n expansion agree where they 
overlap: eq. (7.33) to order (n^^, /c^) is the same as the expansion of the result of [9, 21] 
to the same order. 

Eq. (7.33), however, also contains terms oc k'^ , k^ which are new. They should give 
a more accurate dispersion relation at values /c ~ 1, where the hydrodynamic methods 
become less precise. This accuracy is apparent in figure 1. While eq. (7.33) gives a rather 
poor approximation to the dispersion curves for n = 1,2 (which is not too surprising), 
it gives a very good match already for n = 3. Overall, eq. (7.33) is a much better fit 
to the entire curves than the hydrodynamic approximation in [21]. E.g., for n = 3 the 
hydrodynamic curve (which is exact in n) is better only at relatively small k (< 0.25). 

The agreement between (7.33) and the numerical data is, for large portions of the 
curves, quite better than the expected error ~ 1/n^. The largest deviations tend to appear 
near the threshold mode: for n = 3 we find y /^(^^'^'y**'^) ~ 0.96, which is within the 

~ margins. At n = 4 the match is indeed much better, jS^^™^) jj^^^^-^v^'^^) ~ 0.99. 



^See appendix C for some additional comparison. 
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Figure 1: Dispersion relation of unstable modes for n = 3, 4, 5, 6, 7: the solid line is 
our analytic approximation eq. (7.33); the dots are the numerical solution (the same as in 
[9] , courtesy of P. Figueras) . 

The accuracy at much larger values of n is at least as good as the precision of the 
numerical data we have (~ 10~^ — 10~^). This is shown in fig. 2, where we focus on the 
region near the zero mode where the discrepancies are possibly largest. For reference we 
also include the analytical result from [21]. 

7.5 Final remarks 

AdS black branes at large D. Ref. [21] describes a map that relates the dynamics of 
vacuum black p- branes inD = n + p + 3 dimensions to the dynamics of AdS black branes 
m d + 1 dimensions, by taking n o —d. So, by a suitable analytic continuation of n into 
the region of large negative numbers, we get a map of the solution in this section into a 
solution for perturbations of AdS^+i black branes at large d. Under this map the unstable 
mode becomes a stable, damped, quasinormal mode. 

This solution goes beyond the hydrodynamic regime by being exact in the wavenumber 
k = k/y/d at each order in 1/d. Thus, it contains information about transport coefficients 
at all hydrodynamic orders. In particular, since to the order l/d^ that we have solved, the 
terms fc^ and k^ are non-zero, this gives non-zero values for certain combinations of third 
and fourth order transport terms. 
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Figure 2: Dispersion relation Q{k) for n = 50, 100 for k near the threshold zero-mode. The 
solid line is our analytic approximation eq. (7.33); the dots are the numerical solution; the 
dashed line is the analytic approximation of [21]. 

Scalings. The modes relevant in this section are localized in the sphere of influence at 
r — ro < ro/D. However, the scaling of their frequencies and wavenumbers 

17~0(Z)°), k^O{D^/'^), (7.35) 

is different than the scaling u ~ 0{D) of the characteristic black hole frequencies discussed 
in sees. 4 and 6. They are also the scalings of a specific set of quasinormal modes of AdS 
black branes by the mapping above. Moreover note that k ~ see (2.13), is a scale 
that has not played any role in sees. 4 and 6. 

We suspect that the explanation for this scaling of black brane modes is related to the 
fact that they have a hydrodynamic limit. That is, these are fluctuations along a horizon 
that extends infinitely in one direction, while in our discussion, particularly in sec. 2.3, we 
made a point of taking the limit D — t- oo while keeping the horizon lengths finite. This 
draws attention to the fact that when there are widely separate scales along the horizon 
(at fixed D), the dynamics at both large D and large separation of horizon scales can 
depend significantly on the order in which limits are taken. This issue deserves further 
exploration. 

To conclude, the results in this section are quite encouraging for applications of the 
large D expansion — both because the calculations can be carried out explicitly up to a 
fairly high order, and because of their excellent quantitative agreeement down to rather 
low dimensions. 
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8 Conclusion 



Any gravitational problem that can be formulated in an arbitrary number of dimensions is 
susceptible to study in a large D expansion. Since some terms in the equations will drop 
out of the expansion, this is virtually guaranteed to yield some simplification. Problems 
that otherwise require numerical study may in this approach become analytically tractable. 

The basic physical intuition for why and how this can be so, lies in the remarkably 
consistent picture we have found for General Relativity at large D. To sum it up: 

• In the limit D — )• oo, the spacetimes are flat Minkoswki geometries with 'holes' 
removed at the location of singular horizons. The radii of these holes are finite, but 
the of their horizons vanishes. 

• When — )• oo, the holes do not interact with each other, nor with waves of any 
finite frequency — the horizons are reflective walls of infinitely high potential. Only 
for frequencies and wavenumbers that diverge with D can the interactions remain 
finite in the limit D — )• oo. 

• At large but finite D the black hole has non-trivial dynamics concentrated within a 
thin region that extends out to a distance of order 1 / D from the horizon. 

• This dynamics can be reproduced by an effective theory of fields propagating in flat 
spacetime, subject to specific boundary conditions on the horizons, like (6.56) for 
massless scalar fields. 

We have suggested an image of black holes as dust particles as capturing some features 
of this large D limit, in particular the absence of interactions, the vanishing cross sections, 
and the limiting equation of state of black branes. A different image of large-D black holes 
as fluid lumps was suggested in [8]. It retains the property that the radius of the black 
holes is finite, and that when small 1/D corrections are included, black strings (~ fluid 
tubes) have a Plateau-like instability and break apart into black holes (~ droplets); black 
holes only interact, and then can merge, when their surfaces are within very close distance 
of each other (like fluid drops do). It also suggests an analogy between droplet and black 
hole evaporation when molecular /quantum mechanical effects are included. Although we 
do not have an accurate and complete mapping of the large n black hole physics to either 
of these systems, the images still may, within their limitations, be of some intuitive help. 

One may, perhaps, feel uneasy about regarding a quantity as basic as the number of 
dimensions as a tunable parameter of a theory. However, this is not different than in 
SU{N) gauge theories, or actually in most theories with adjustable parameters. We are 
viewing the gravity theories at different values of D as being one uniparametric theory, 
and eq. (2.1) as one uniparametric solution of this theory. 

We have studied only one of the limits of the parameter space of the theory, D — t- oo, 
but it may also be interesting to study the opposite limit. Arguably this corresponds not 
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to = or 1 but to D = 3. At D = 3, the theory (1.1) ceases to have both local dynamics 
and black hole solutions, but taking D as a continuous parameter one still finds non-trivial 
effects for D = 3 + e. This idea was put forward in ref. [7] to solve the problem of the 
Schwarzschild negative mode in an expansion in e. It should be of interest to understand 
this expansion in more generality. ^'^ 

Our study of gravitational radiation has been rather cursory, since at — )• oo it 
largely decouples from black holes — except in phenomena where the radiation probes 
scales extremely close to the horizon. Nevertheless, gravitational waves are a sector of the 
theory that deserves closer attention. Some features of classical gravitational radiation 
depend on the spacetime dimensionality, e.g., on whether D is even or odd, as is the case 
for the applicability of Huygens principle [48, 39] . Further work is required to clarify the 
importance of all such effects in the large D expansion. 

We have not considered the coupling of the gravitational field to other fields or matter 
systems, except for brief references to charge and the cosmological constant. The main 
issue when considering these additions is the choice of how the new parameters scale with 
D. What choice is most convenient will depend on the particular type of problem one 
is interested in, but often it is useful to keep fixed the length scale that characterizes 
the effects of the new fields on the geometry. These remarks apply as well to solutions 
with compactified dimensions. The size of the compact space may be fixed/grow/shrink 
as D increases, each choice capturing different physics. In a similar vein, in the context 
of higher-dimensional gravity it is natural to also include the class of Lovelock theories. 
Each of the new terms that these theories introduce comes with a new length parameter. 
Depending on how these are chosen to scale with D, different large D limits will result. The 
fact that the Riemann curvature tends to strongly localize close to the horizon indicates 
that the dust mote picture should still apply at least in some situations {e.g., the 'normal 
branch' of Lovelock black holes). 

Finally, another setting in which the large D expansion might be useful is the study of 
cosmologies of four-dimensional braneworlds in a large D bulk. Having fixed worldvolume 
dimension, the gravitational effect of the brane in the bulk will be strongly localized close 
to the brane, which may simplify some calculations. 
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A Elementary geometry at large D 

Consider D-dimensional Minkowski spacetime. Using 

/ de{si^ef = ^ , (A.l) 

one obtains the area of the unit-radius sphere 5^~^ 

Aspi,= [ dOi [ d92... [ dOo-z I #(sin^i)^-3(sin02)^-^..sin0D-3 
Jq Jo Jo jo 

= (A.2) 

as given in (2.3). This is not a monotonic function of D: it reaches a maximum at 
D ~ 8.257, and then rapidly decreases as D grows. For integer D the maximum is at 
D = 8, where = 167r'^/15 ~ 33. One might be tempted to conclude that this non- 
monotonicity could imply qualitative differences between values of D smaller or larger 
than 8, and thus a potential inadequacy of the large D expansion when applied to D < 8. 

However, this result does not make clear what we are comparing the sphere to. A more 
careful analysis shows that in an appropriate sense spheres do become monotonically small 
as D grows. The area of the circumscribed cube that contains the unit-radius sphere and 
is tangent to it at its faces, is 

^circubc = {D- 1)2^-1 . (A.3) 

The ratio AsphMcircube can now be seen to be a monotonically decreasing function of 
D > 2, which vanishes as 



^circube ^ 2Z) / 



0. (A.4) 



A simple intuition for this behavior follows from considering that the length L of the 
diagonal of a cuboid in D spacetime dimensions is 

= X? + • • • + . (A.5) 

Then, for a generic cuboid at large D, its diagonal length L is much larger than any of its 
side lengths Xj. In particular, the length L of the half-diagonal of the circumscribed cube 
around a unit-radius sphere is 

L = - 1 > 1 . (A.6) 
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Now a plain indication that the sphere is in fact semifactorially smaller than its circum- 
scribing cube is that 

^ ^ D-""^ , (A.7) 



L 

which reproduces the leading behavior in (A. 4). 

If we consider the inscribed cube (with its vertices on the sphere) then 

Ancube = {D- 1)2-^/22^-1 , (A.8) 

which becomes much smaller than the sphere 

^^^J^y^^O. (A.9) 

The behavior is monotonic for D > 2, but the shrinking rate is slower than in (A. 4), since 
Ancube ~ D~^/'^, as follows from the above argument about the diagonals. 

Similar behavior is obtained by comparing the volumes V inside these bodies, since 
for the unit-radius sphere and for the cube of side-length two we have 

V = ^. {A.10) 

However, these volumes do not appear to be relevant for any purpose in this paper. 



B Technical appendices to sections 6 and 7 

B.l Far scalar wave solution in the overlap region 

In the overlap region we have In R ^ n and so we can expand 

r = l + -lnR + 0(n-^). (B.l) 

n 

Since at large n and fixed r both the argument and the index of the Bessel functions 
in (6.34) are large numbers of the same order ~ n,^^ it is appropriate to use the Debye 
expansion. When uj < Uc it gives 

V'(r) = ^ = /'cie-"'^^("~*^'"^") - 2C2e"^^("-*^''^") +0(n~i)) , (B.2) 

^y2^TRnuJctanha V ^ 

while at cD > Wc, 



i'{r) = ^ \ ((C7i-iC2)e--""-(^-^'^"^)---/^ 



+ {Ci + iC2)e^"^^(^-*^°'^)+*^/^ + O(n-i) ) . (B.3) 



*So, again, the expansion is not valid when lj > n. 
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Here a and (3 are defined by 



— r 



sech a , lu < Uc , 
sec /3 , < a) . 



(B.4) 



Expanding in tlie overlap region using (B.l) we find 

InR 



a — tanh Q = Qo — tanh ao — tanh ao \- 0(n ), (B.5) 

n 



and 



In R 

/3 - tan/3 = /3o - tan/3o - tan/3o h 0(n"^) . (B.6) 

n 

With this, the far solution in this region can be written as in eq. (6.37). 
B.2 Black brane far solution at next order 

In the far region r~" is exponentially small in n and therefore it does not yield any 
perturbative 1/n corrections. In order to compute corrections to the leading far solution 
(7.9) we introduce an auxiliary order-counting parameter e, 

fir) = 1 - er-" , (B.7) 

which is set to e = 1 at the end of the calculations. The far solution is then expanded in 
the form 

The first order equation in the far region is 



ar^ r ar \ r'^ J \ 



^^^^^°Vg(V°M, (B.9) 



dr 



where 



Green's method gives the solution as 



(1) ^ ^a)^^^ M^ + r,^^, (B.ii) 

' j.n/2 j.n/2 j.n/2 ' ^ ' 

where 

Si= K^{x)S{x)dx, ^1=/ Iuix)Six)dx. (B.12) 
5 is a source term defined as 

5 = -fe^V^ f ^^'^^ + Q'^'^'^^A ■ (B.13) 
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We can take f arbitrarily since it does not affect the matching procedure. 

A very lengthy calculation yields Si and Ti in an unilluminating form. For the purposes 
of matching one only needs the expressions in the overlap region. With the normalization 
of the leading order solution (7.9) fixed like in sec. 7.2 (eq. (7.12)), namely 

then the solution in this region is 

+^ — ( Pn^ -n^ + k^ + 2k^n^ + 3p - k^ 



2A;4n2R2 

+ {k^ -Q'^- SPn"^ + 2k^ + P) In R + k^{l + A:^)(ln R)^ ' 
+0(n-3). (B.15) 

The actual expression for An is rather complicated (dependent on n, k and 17) and we 
omit it. 

This solution yields terms oc R^^ ^t orders n^, n^^ and n^'^. We can match them to 
the corresponding terms in the large R expansion of the near solutions r/(o) (7.24), r/(i) 
(7.25), and 7/(2) (7.28). The first of these is just a matching of the overall amplitude. The 
second one and the third one would provide, respectively, the dispersion relations (7.27) 
and (7.32). 

B.3 Sources for near-region equation 

The source term of the n-th order equation is given by 

n 

where 7?(j) is the j-th order solution. The S'(j) at each order are 

5(1) (r/) = -2(R - if [1 + 2R(-1 + 2R)P] r?' 

- (R- 1)2[3 + (1 -4R + 8R2)P]r/, (B.17) 



2 



5(2)(r/) = -8(R - 1)3r[RJ72 + (_i + r + (_i + 2R) In R)^ 

+ 2R(l-3R + 2R2)P]r/' 
- [1 - 2R + R2(1 - 5n^) - 402r3 + 8fi2R4 

+ 2R(R - 1)(4 - 12R + 8R2 + (1 - 4R + SR^) In R)P 

+ 8R2(R- 1)2(1 - 2R + 4R2)A;%, (B.18) 
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and 



5(3)(?7) = -8(R-1)^R J7^R(1 + 21nR) + r(2R^(-3 + 4R)fi^ 
+ 2(-l + R) In R + (-1 + 2R)(ln R)^ 
+ 4R(R - 1)P(-1 + R + (-2 + 4R) In R + 2R(R - 1)(-1 + 2R)P))j ry' 

- 2(R - 1) [Rr?2(2(R - l)(-l + 4R) + (-5 - 4R + SR^) In R) 
+ (R - 1)P(4(1 + R(-2 + R + R(-l - 4R + SR^)^^^)) 
+ 8(R - 1)(-1 + 2R) In R + (1 - 4R + 8R2)(ln R)^ 
+ 16R(R - 1)P(1 - 3R + 2R2 + (1 - 2R + 4R2) In R 

+ R(R-l)(l-2R + 4R2)P))jr/ (B.19) 

5(4)(r?) = -^^^^^^^^[3RJ^2(2R2j^2^1nR + (lnR)2) + P(12172R2(R_ 1) 
3 

+ 12R2(-3 + AR)n^ In R + 3(R - l)(ln R)^ + (-1 + 2R)(ln R)^ 
+ 24(R - 1)RP(R2(-2 + 3R)n^ + (-1 + R + (-1 + 2R) In R) In R 
+ (R - 1)RP(-1 + R + (-3 + 6R) In R + 2(R - 1)R(-1 + 2R)P)))j rj' 

r3Rj^2(4^p _!)(_! + R + 2R2(1 + 2R)n^) + 4(R - 1)(-1 + 4R) In R 

3 

+ (-5 - 4R + 8R2)(ln R)^) + 2(R - 1)P(12(R - 1)R2(-3 + 8R)n^ 
+ 12(1 + R(-2 + R + 2R(-1 - 4R + 8R'^)n'^) In R 
+ 12(R - 1)(-1 + 2R)(ln R)2 + (l _ 4R + 8R2)(ln R)=^ 
+ 24(R - 1)RP((R - 1)2 + 6R3(-1 + 2R)n^ 
+ 2 In R(2 - 6R + 4R^ + (1 - 2R + 4R2) In R) 
+ 2(R - 1)RP(2 - 6R + 4R2 + 3(1 - 2R + 4R2) In R 
+ 2(R - 1)R(1 - 2R + 4R2)P)))j rj. (B.20) 

Sources at large R. The integrations constants Ai and Bi at i-th order in (7.17) are 
determined by the boundary conditions at R = 1 and R S> 1. This i-th order solution 
contributes to the large R condition for the {i + l)-th order solution through ^(^^(//(j)). 
Since (rj) is linear in r], the contributions of Ai and Bi can be obtained independently. 
Now, 

S(^i){AiUo + Bivo) = -Ai{l + P)R - Bi(8R2 - 12R)P + 0(R°), (B.21) 

so the leading term at large R is controlled by Bi, not Ai. This source term is contained 
in the condition at large R as 

1 

^ no(R')S(i)(Aino + 5it;o)(iR' = -4p5i + 0(R-i), (B.22) 

i.e., only Bi enters in the large R condition for the [i + l)-th order solution. This is why 
in sec. 7.3 we can proceed from the third to the fourth order without A^. 
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C Hydrodynamic vs. large D expansions 



When refs. [9, 21] compared the hydrodynamical and numerical calculations of Q{k), it 
was observed that the agreement improves for larger n. The proposed explanation for 
this relies on a conjecture about the large n-dependence of higher-order hydrodynamic 
transport coefficients. More precisely, the large D expansion in this article is of the form 

n = + ^ + (C.I) 

In the hydrodynamic expansion it seems more appropriate to keep fixed T = n/(47rro), 
instead of rg. Rescaling VL = {n + k = \/n + 1 A:, the expansion is 

"-'('-^4lr + -("'(4lf)' + -(»'(4lf)' + -) ■ (") 

The expansion variables in the two cases are equivalent at large n, since then k = k/{4:TrT). 
Nevertheless, these are clearly different expansions. It is apparent that the large-n accuracy 
of the hydrodynamic results requires that Cj{n) — )• when n — >• oo for all j [9]. Hence, 
lacking an explicit calculation of, say, C2(n), we could not be sure within the hydrodynamic 
approach (i.e., prior to comparing to the numerical results) that at large n there will not 
remain a term oc k^ in the dispersion relation at leading order. 

The Cj{n) are obtained from effective hydrodynamic transport coefficients computed 
from black brane perturbations. While it seems possible that their leading large n scal- 
ing can be determined from generic features of higher-dimensional gravity, the required 
behavior has not been derived yet within the hydrodynamic approach. Our calculation of 
the dispersion relation in the large n expansion can be regarded as a proof of it. 

The hydrodynamic expansion applies also for non-linear black brane perturbations. It 
should be interesting to investigate if the large D expansion can be useful in that problem 
too. 
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